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We address conditional de-Gaussification of continuous variable states by inconclusive photon subtraction 
(IPS) and review in details its application to bipartite twin-beam state of radiation. The IPS map in the Fock 
basis has been derived, as well as its counterpart in the phase-space. Teleportation assisted by IPS states is 
analyzed and the corresponding fidelity evaluated as a function of the involved parameters. Nonlocality of IPS 
states is investigated by means of different tests including displaced parity, homodyne detection, pseudospin, 
and displaced on/off photodetection. Dissipation and thermal noise are taken into account, as well as non unit 
quantum efficiency in the detection stage. We show that the IPS process, for a suitable choice of the involved 
parameters, improves teleportation fidelity and enhances nonlocal properties. 



I. INTRODUCTION 

Nonclassical properties of the radiation field play a rele- 
vant role in modern information processing since, in general, 
improve continuous variable OC V) communication protocols 
based on light manipulation jll 12]. Indeed, quantum light 
finds application in several fundamental tests of quantum me- 
chanics 1 3], as well as in high precision measurements and 
high capacity communication channels |5|. Among non- 
classical features, entanglement plays a major role, being the 
essential resource for quantum computing, teleportation, and 
cryptographic protocols. Recently, CV entanglement has been 
proved as a valuable tool also for improving optical resolution, 
spectroscopy, interferometry, tomography, and discrimination 
of quantum operations. Recent experimental realizations also 
include dense coding |6] and teleportation network |7j]. 

Entanglement in optical systems is usually generated 
through parametric downconversion in nonlinear crystals. The 
resulting bipartite state, the so-called twin-beam state of radi- 
ation (TWB), allows the realization of several beautiful ex- 
periments and the demonstration of the above quantum pro- 
tocols. However, the resources available to generate CV en- 
tangled states are unavoidably limited: nonlinearities are gen- 
erally small, and, in turn, the resulting states have a limited 
amount of entanglement and energy. In this context, practical 
applications require novel schemes to create more entangled 
states or to increase the degree of entanglement of a given sig- 
nal. 

In quantum mechanics, the reduction postulate provides an 
alternative mechanism to achieve effective nonlinear dynam- 
ics. In fact, if a measurement is performed on a portion of a 
composite system the output state strongly depends on the re- 
sults of the measurement. As a consequence, the conditional 
state of the unmeasured part, i.e. the sub-ensemble corre- 
sponding to a given outcome, may be connected to the initial 
one by a (strongly) nonlinear map. In this paper, we focus our 
attention on a scheme of this kind, and address a conditional 
method based on subtraction of photons to enhance nonclas- 
sical features. In particular, we analyze how, and to which 
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extent, photon subtraction may be used to increase nonlocal 
correlations of twin-beams. As we will see, photon subtrac- 
tion transforms the Gaussian Wigner function of TWB into 
non-Gaussian one, and therefore it is also refeiTed to as a de- 
Gaussification process. 

The photon subtraction process on TWBs was first pro- 
posed in where a well defined number of photons is be- 
ing subtracted from both the parties of a TWB, by transmit- 
ting each mode through beam splitter and performing a joint 
photon-number measurement on the reflected beams. The de- 
gree of entanglement is then increased and the the fidelity of 
the CV teleportation assisted by such photon subtracted state 
is improved |9]. However, this scheme is based on the possi- 
bility of resolving the actual number of revealed photons. In 
1 10] we showed that the improvement of teleportation fidelity 
is possible also when the number of detected photon is not 
known. In our scheme we use on/off avalanche photodetec- 
tors able only to distinguish the presence from the absence of 
radiation. For this reason we referred to this method as to in- 
conclusive photon subtraction (IPS). The single-mode version 
of this process has been recently implemented LI lil and the 
nonclassicality of the generated state starting from squeezed 
vacuum has been theoretically investigated | 3 fl3ll . In ad- 
dition, nonlocal properties of the photon-subtracted TWBs 
have been investigated by means of different nonlocality tests 
15, 16, 17, 18, 19], finding enhanced nonlocal proper- 
ties depending on the particular test and on the choice of the 
involved parameters. 

This paper is devoted to review the effects of IPS process on 
TWBs either in the ideal case, i.e., when the detection are not 
affected by losses and no dissipation or thermal noise occurs 
during the propagation of the involved modes, or when non 
unit quantum efficiency is taken into account as well as the 
dynamics through a noisy channel is considered. 

The paper is structured as follows: in the next section we 
introduce photon subtraction as a method to enhance nonclas- 
sicality of a radiation state and illustrate inconclusive photon 
subtraction on a single-mode field. The de-Gaussification pro- 
cess on two-mode fields is described in Sec. |llll where the 
map of the IPS process is given both in the Fock representa- 
tion and in the phase-space. In Sec. lIVI we briefly review the 
dynamics of a TWB in noisy channels and show that IPS can 
be profitably applied also in the presence of noise. The CV 
teleportation protocol is described in Sec.|VJ where we com- 
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pare the teleportation fidelity when the protocol is assisted or 
not by the IPS process. In the following Sections, in order 
to characterize in details the nonlocal properties of the IPS 
states, we consider different Bell tests, namely, the nonlocal- 
ity test in the phase space (Sec. lVH . the homodyne detection 
test (Sec. lVlfl . the pseudospin test (Sec. IVIlTl l. and a nonlo- 
cality test based on on/off photodetection (Sec. lIXI i. Finally, 
Sec. Incloses the paper with some concluding remarks. 



II. PHOTON SUBTRACTION 
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FIG. 1: Scheme of the IPS process: the input state Qs is mixed with 
the vacuum state go = |0){0| at a beam splitter (BS) with trans- 
missivity T; then, on/off photodetection with quantum efficiency e is 
performed on the reflected beam. When the detector cUcks we obtain 
the IPS state gi 



The idea of enhancing nonclassical properties of radiation 
by subtraction of photons has been introduced in the context 
of Schrodinger cat generation |20] and subsequently applied 
to the improvement of CV teleportation fidelity |8]. In the 
schemes of Refs. |8, 20] the field-mode to be "photon sub- 
tracted" (PS) is impinged onto a beam-splitter with high trans- 
missivity and whose second port is left unexcited. At the out- 
put of the beam splitter the reflected mode undergoes pho- 
ton number measurement whereas the conditional state of the 
transmitted mode represents the PS state. The properties of 
the PS state depend on the number of detected photons, with 
single-photon subtracted states that play a major role in the 
enhancement of nonclassicality. Unfortunately, the realiza- 
tion of photon number resolving detectors is still experimen- 
tally challenging, and therefore a question arises concerning 
the experimental feasibility of subtraction schemes. 

Photodetectors that are usually available in quantum optics 
such as avalanche photodiodes (APDs) operates in the Geiger 
mode 1 21, 2 2) . They can be used to reconstruct the pho- 
ton statistics 11231 l24ll but cannot be used as photon counters. 
APDs show high quantum efficiency but their breakdown cur- 
rent is independent of the number of detected photons, which 
in turn cannot be determined. The outcome of these APD's is 
either "off" (no photons detected) or "on", i.e., a "click", in- 
dicating the detection of one or more photons. Actually, such 
an outcome can be provided by any photodetector (photomul- 
tiplier, hybrid photodetector, cryogenic thermal detector) for 
which the charge contained in dark pulses is definitely below 
that of the output current pulses corresponding to the detection 
of at least one photon. Note that for most high-gain photmulti- 
pUers the anodic pulses corresponding to no photons detected 
can be easily discriminated by a threshold from those corre- 
sponding to the detection of one or more photons. 

It appears therefore of interest to investigate the properties 
of photon subtracted states when the number of detected pho- 
tons is not discriminated. Such a process will be referred to 
as inconclusive photon subtraction (IPS) throughout the pa- 
per. The scheme of the IPS process is sketched in figure ^ 
The mode a, excited in the state Qs is mixed with the vacuum 
qq = |0)(0| (mode b) at an unbalanced beam splitter (BS) 
with transmissivity T = cos^ <j) then, on/off avalanche 
photodetection with quantum efficiency e is performed on the 
reflected beam. APDs can only discriminate the presence of 
radiation from the vacuum. The positive operator-valued mea- 



sure (POVM) {no(e), Hi (e)} of the detector is given by 



no(e) 



fe=0 



{l-e)''\k)(k\, ni(e)=I-no(e). (1) 



The whole process can be characterized by T and e which will 
be referred to as the IPS transmissivity and the IPS quantum 
efficiency. The conditional state of the transmitted mode after 
the observation of a click is given by 
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-Tr. 



Uabi(f>)gs «) go C/]b(0) la ® ni(e) 



(2) 



where Uab{4>) — exp{—(j){a''b — a6^)} is the evolution oper- 
ator of the beam splitter, and pi{(f>, e) is the probability of a 
click. In general, the transformation (|2j realizes a non unitary 
quantum operation gi = £{gs) with operator-sum decompo- 
sition given by 
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J2mpicly,e)Ep{cP)g,El{(b) (3) 



where 



tan^P [1 - (1 - s)P] 



Mp{(t)) = flP cos" "^0. 



(4) 
(5) 



which is found by explicit evaluation of the partial trace in (|2j. 
The IPS state obtained by applying the map (|3} to a Gaussian 
state is no longer Gaussian, and therefore IPS represents an 
effective source of non Gaussian states, which should be oth- 
erwise generated by highly nonlinear, and thus inherently low 
rate, optical processes. 

In general the IPS process can produce an output state 
whose energy is larger than the one of the input state and 
whose nonclassical properties can be enhanced. As an ex- 
ample, we address the photon subtraction onto a Gaus- 
sian state described by the following Wigner function (using 
the Wigner function formalism makes analytical calculations 
more straightforward): 



exp{— 



G{z^ + z*^)} 



(6) 
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whose energy is given by 



F 1 
(7) 



When the state (|6j undergoes the IPS process described above, 
the Wigner function associated with the output state qi reads 



k=l 



with Ci{(j),e) = 1, C2(0,e) 
where 



B = (l-T)cr+|l2, (TM 



^ ^C,(0,e)W^fc(z), (8) 
2-e 



2£ 



l2 



(9) 



I2 being the 2x2 identity matrix, and cr is covariance asso- 
ciated with the state (|6|l 



{F + 2Gy^ 

(F-2G)-i 



(10) 



where [aU = ^{{Rh,Rk}) - {Rh){Rk), {A, B} = AB 
BA denotes the anticommutator, and 



R — {Ri, R2) 



a + flT a — a' 



\/2 ' 



(11) 



(• • • being the transposition operation. Notice that Wi{z) 
is no longer Gaussian. In Eq. (|8j we defined 



Wk{z) 



where 
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Gi = -(Z^+-Z^_), 



with 



V± 



F2 = 2(V+ + V_) , G2 = V+ - V_ , 
F±2G 



2T+ (1 -T)(F±2G) ' 

F + 2(1 ±G)T 
4r + (1 - T)[2e + (2 - £)(F ± 2G)]} 



(12) 

(13) 
(14) 

(15) 
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Because of the analytical expression (|8), the energy of the 
photon subtracted state is simply given by 



Ei{T,e) 
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FIG. 2: Logarithmic plots of the energies (dashed line) and Ei 
(solid lines) in the case of a squeezed vacuum |0, r) as input state 
and as functions of tanh r for e = 1 and different values of T. From 
top to bottom (solid lines): T = 1, 0.9, 0.75, and 0.5. 




with Ck = Ck {T, e) and where we put T = cos^ <j). 

Let us now focus our attention on the IPS process applied to 
the squeezed vacuum 1 0, ?■) = S{r)\0), S{r) = exp{ir(a^ — 
a^)} being the squeezing operator, which has been recently with equal transmissivity T — cos^ (p; the reflected modes 
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FIG. 3: Plots of the energy Ei of the IPS state in the case of a 
squeezed vacuum j0,r) as input state as a function of tanhr for 
T = 0.9 (solid lines) and T — 0.8 (dashed lines) and different val- 
ues of e. From top to bottom: e = 1, 0.75, 0.5, and 0.25. 



realized experimentally fT]\. The Wigner function associated 
with |0, r) is given by Eq. (|6j with F ~ 2cosh2r and G ~ 
— sinh2r. In Figs. |5]we plot the energies Eg and Ei of the 
input and output states, respectively, for different values of 
the involved parameters as functions of tanh r. We can see 
that there is a threshold on r, depending on T and e, under 
which the IPS state has a larger energy than the input state. 
Furthermore, when e = 1, T — ^ 1 and r — > we can see 
that El — > 1: in these limits the output state approaches to 
the squeezed Fock state S'(r)ll) L12., .131 . Finally, in Fig.|3 
El is plotted for two values of T and different values of e 
as a function of tanh r: we find that as r increases, the IPS 
efficiency is not so relevant in the process. 



III. PHOTON SUBTRACTION ON BIPARTITE STATES 

In this Section we address de-Gaussification of bipartite 
states by IPS . The de-Gaussification can be achieved by sub- 
tracting photons from both modes through on/off detection 
llgl fioll . The IPS scheme for two modes is sketched in Fig.0] 
The modes a and b of the shared bipartite state Qs are mixed 
with vacuum modes at two unbalanced beam splitters (BS) 
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FIG. 4: Scheme of the IPS process. The two modes, a and b, of a 
shared bipartite state Qs are mixed with the vacuum at two BSs with 
equal transmissivity T and on/off photodetection with quantum effi- 
ciency e is performed on the reflected beam: when both the detectors 
click one obtains the IPS state. 



J 



and d are then revealed by avalanche photodetectors (APD) 
with equal efficiency e. The conditional measurement on 
modes c and d, is described by the POVM (assuming equal 
quantum efficiency for the photodetectors) 



noo(e) = no,c (e)(8) no^d(e), 
noi(e) = no,c(e)(8)ni^d(e) , 

nio(e) = ni,e(e)(8)no^d(e), 
nii(e) = ni,e(e)(8)ni,d(e) . 



(18) 
(19) 
(20) 
(21) 



When the two photodetectors jointly click, the conditioned 
output state of modes a and b is given by M\1M 



^, . Tr,rf [Uac{4') <» Um{4') Qs <8 \0)cddc{0\ Ul{(j)) (8 la^h^ 

Pii{i;(P,e) 



nii(e)] 



(22) 



where Uac{(f>) = exp{- 



1^)} and [/fcc;(0) are the evo- with 



lution operators of the beam splitters, \0)cd = |0)c <8i |0)d, and 
Pii(r, (j>, e) is the probability of a click in both the detectors. 
The partial trace on modes c and d can be explicitly evaluated, 
thus arriving at the following decomposition of the IPS map: 



X J2 rnp{cp,e)Mpq{4>) Q,Ml^{4>)mq{4>,e) (23) 

p,q=l 



where 



(24) 



Eq. ( I23> is indeed an operator-sum representation of the 
IPS map: {p, q} = 9 should be intended as a polyin- 
dex so that i23\ reads £{gs) = J^o ^eQsA^ with Ag — 
[pii(r, s)]-^/^ mp{(l), e) Afp,(0). 

From now on we focus our attention on the case in which 
the shared state is the twin-beam state of radiation (TWB) 
= |A))((A|, where |A)) = ^/T^J2^ \''\k) (E) \k) with 
A — tanh r, r being the TWB squeezing parameter The TWB 
is obtained by parametric down-conversion of the vacuum, 
|A)) = exp{r(a^&^ — ab)}\0), a and b being field operators, 
and it is described by the Gaussian Wigner function 



Wo(a,/3) 



exp{-2Ao(|a|' 



+ 2Bo{af3 + a* (3*)} 



^VDet[o-o] 



(25) 



An 



4v'Det[o-o] 



So 



4v/Det[cro] ' 



(26) 



where Aq = Ao{r) = cosh(2r), Bq = Bo{r) = sinh(2r) 
and (To is the covariance matrix 



era = 



1 I A0I2 Bo(T3 

2 I Bo 0-3 Ao 1-2 



(27) 



I2 being the 2x2 identity matrix and 0-3 = Diag(l, —1); (Tq 
is defined as [cro]hk = ^({Rh,Rk}) - {Rh){Rk) with 



-R — (i?i, i?2, RsiRi)^ 

'a + a"f a -at b + b^ b - b'^ 
a/2 ' zV2 ' ^/2 ' i\/2 



(28) 
(29) 



Now we explicitly calculate the Wigner function of the corre- 
sponding IPS state, which, as one may expect, is no longer 
Gaussian and positive-definite. The state entering the two 
beam splitters is described by the Wigner function 

Tn'°'(«,/?,C,a = W^o(«,/3)^exp{-2|Cp-2|^|2} , 

(30) 

where the second factor at the right hand side represents the 
two vacuum states of modes c and d. The action of the beam 
splitters on W,'."" can be summarized by the following change 
of variables |2] 



a a cos 4> + C sin cf) , 
P ^ f3 cos ^ + ^ sin (/) , 



^ ^ ^cos^ — asin0, (31) 
^ ^ £, cos (p — PsiiKp , (32) 
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and the output state, after the beam splitters, is then given by 

w^:;'(«,/3,c,e) = 

^WrAa,f3) exp{-a|e|'+we + w*r} 
X exp { - a|Cp + iv + 2Bo^ sin^ 0)C 

+ (t;*+2Borsin2 0)C}, (33) 

where 

exp{-6(|a|2 + |^|2) + 2Bo cos^ ^ (a/? + a* (3*)} 



7r2\/Det[cro] 



(34) 



and 



a = a(r, 0) = 2(Aosin^(/) + cos^(?!)), (35) 

6 EE 6(r, (/)) = 2(lo cos^ + sin^ 0) , (36) 

V = v{r,(t)) = 2cos(/) sin (/>[(! - Ao)a* + BqP], (37) 

u; = u;(r, 0) = 2 cos </) sin [(1 - Ao)/?* + Boa] • (38) 

At this stage on/off detection is performed on modes c and 
d (see Fig. |4}. We are interested in the situation when both 
the detectors cHck. The Wigner function of the double click 
element nii(e) of the POVM [see Eq. (EB] is given by uH 

m 



(39) 



^{1 - ge(C) - QeiO + QeiOQeiO} , (40) 



with 



Qeiz) 



2-e 



exp ■ 



2e 

2-e 



(41) 



Using Eq. i22i and the phase-space expression of trace for 
each mode, i.e., 

Tr[Oi 02]^7r f d^z W[Oi]{z) W[02] (z) , (42) 
Jc 

the Wigner function of the output state, conditioned to the 
double click event, reads 



/(a,/3) 



pii(r,^,e) ' 
where f{a,P) = /r,0,e(a,/3) with 



(43) 



C2 TT" 

4 



E5G?ie(«-/3-C,e), (44) 



k=l 



with Ck = Ck{e) and Ci = 1, C2 = C3 = -2(2 - e)-\ 
C4 = 4(2 — e)^^; the double-click probability pn (r, 0, e) 
can be written as function of /(a, /?) as follows 



pii(r,0,e) =7r^ / d'a d' P f{a, (3) . (45) 
The quantities G^^l ^(a, /3, C, in Eq. are given by 



exp 

+ {v* + 2BC sin^ 0)C*} 

xexp{-yfc|eP+< + w;*r} , (46) 

where Xk = Xk{r, (j), e), and yt = yk{r, </>, e) ai-e 

a;i = = yi = y2 = a 

X2 ^ Xi ^ y^i = Vi ^ a + 2e{2 - e)"^ . 

After the integrations we have 

1 4 

/(a, /?) = — ^Cfe exp{(/fc - + (5,. - 6)|/3|2 



fe=i 



+ (2BoT + //,fc)(a/3 + a*/3*)} (47) 



and 



pii (r,T,£) = ^ 



t {b ~ fk){h - gk) - {2BoT + hkf 



where we put T = cos'^ (p — 1 ~ siv? (f>, and defined 



Ck = Ck{r, T,e) 



^Ck 



(48) 



(49) 



and fk = fk{r,T), gk = gk{r,T), and hk = hk{r,T) given 
by 

Jk = Mk [xkBl + 4^2(1 _ Io)(i _T)^yk{\- lo)'] , 

(50) 

fffc = A4 [a;fc(l - Mf + 4B2(1 - Io)(l - T) + ^/^S^] , 

(51) 

^fc = Mk {{xk + y/c)5o(l - ^0) 

+ 2Bo[B2 + (1-Io)'](1-T)}, (52) 
4r(l-T) 



Mk^Mk{r,T) = 



Xkyk-^Bl{\-TY 



(53) 



In this way, the Wigner function of the IPS state can be rewrit- 
ten as 

1 4 

W^ips (a, /3) = ^ E ^''^ ^) ' (54) 



7r2pii(r,T,e) ^ 
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FIG. 5: Logarithmic plots of the energies (dashed line) and Eips 
(solid lines) in the case of a TWB as input state as functions of tanh r 
for e = 1 and different values of T. From top to bottom (solid lines): 
T = 1, 0.9, 0.75, and 0.5. 
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FIG. 6: Plots of the energy Eies of the IPS state in the case of the 
TWB as input state as a function of tanhr for T = 0.9 (solid lines) 
and T = 0.8 (dashed lines) and different values of e. From top to 
bottom: e = 1, 0.75, 0.5, and 0.25. 



corresponding energy is 

E = I (fv (fw {\v\^ + |wp - 1) W{v, w) . 



(58) 



If the bipartite state has a Wigner function of the form 

cxp{-F|w|2 - G|wp + H{vw + v*w*)} 



Ws{v,w) 
then its energy reads: 



7r2(FG-iJ2)-i 



F + G 



2{FG-H^ 



1 



(59) 



(60) 



thereby, in the case of a TWB as input state F, G, and H are 
obtained from Eq. i25\ and the energy of the state emerging 
from the IPS process can be written as 



1 



7r2pn(r,T,e) 



k=l 



Fk + Gfe 



2{FkGk - Hl f 



- 1 



(61) 



with Fk = b- fh, Gk = b- gh, and Hk = 2BaT + 
and all the involved quantities are the same as in Eq. ( 154k As 
in the single mode case, we can see that there is a threshold 
on r, depending on T and e, under which the IPS state has 
a larger energy than the input state. In Fig.|6l£'ips is plotted 
for two values of T and different values of e as a function of 
tanh r: we find that as r decreases, the IPS efficiency is not 
so relevant. 

The state given in Eq. \5A\ is no longer a Gaussian state and 
its use in the improvement of continuous variable teleportation 
1101 as well as in the enhancement of the nonlocality 1 14, id 
nM will be investigated in the following Sections. 



with 



IV. DYNAMICS OF TWB IN NOISY CHANNELS 



Wk{a,(3) = exp{(/fe - 6)|a|2 + [g^ - 

+ (2SoT + /ife)(a^ + a*/3*)}. (55) 

Finally, the density matrix corresponding to Wips(a, (5) reads 
as follows fioll 



Mi: 



pii(r, r,£) 



^ (Ar)"+™ A,fc(r,£) 



n,m— 



nnvnra 

X 75;'^7J""^ I" ~ ^/"l"- ~ 'Vbb{m - h\a{m - k\ , (56) 
where A = tanh r and 



AHAT.e)=[l-{l-ef] [l - (1 - e)^ (^-^) 



h+k 



(57) 



In Fig.|5lwe plot the energies Eg and E,ps of the bipartite 
input and output states, respectively, for different values of 
the involved parameters as functions of tanh r. We recall that 
for a given Wigner function W{v,w) of a bipartite state, the 



Before addressing the properties of the IPS bipartite state 
described in the previous Section, we review the evolution of 
the twin-beam state of radiation (TWB) in a noisy environ- 
ment, namely, an environment where dissipation and thermal 
noise take place 1 15]. As we will see, we can include in our 
analysis the effect due to the propagation through this kind of 
channel by a simple change of the involved quantities. Using 
a more compact form, Eq. i25i can also be rewritten as 



WoiX) 



cxp 



X 



(62) 



7rVDet[cro] 

with X = {xi,yi,X2,y2V, a = 75(2^1 + m) and /? = 
-^{^2 + iy2), and (• • • denoting the transposition opera- 



tion. 



When the two modes of the TWB interact with a noisy en- 
vironment, namely in the presence of dissipation and thermal 
noise, the evolution of the Wigner functi on (I25ll is described 
by the following Fokker-Planck equation I25ll26ll27ll 



dtWt{X) 



^ITX 



dlTa^dx]Wt{X), (63) 
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with 9x = {dxi,dy^,dx2,dy2)'^ . The damping matrix is 
given by r = 0^,^;^ rfel2, whereas 



^« 



a- 

where is the 2x2 null matrix and 

I f l + 2Nk 



(2) 



(64) 



2 







l + 2Nk 



(65) 



r^;, A^fe denote the damping rate and the average number of 
thermal photons of the channel k, respectively, cr^ repre- 
sents the covariance matrix of the environment and, in turn, 
the asymptotic covariance matrix of the evolved TWB. Since 
the environment is itself excited in a Gaussian state, the evo- 
lution induced by ( I63l l preserves the Gaussian form ( I62l i. The 
covariance matrix at time t reads as follows O EtIi 



(Tt = G^' (To G;'' + (1 - (Gt) (Too , 



,1/2 



(66) 



where Gt ~ i e l2- The covariance matrix <Tt can 
be also written as 

1 f At{T,,N,)l2 BtiTi)a3 \ 

ft = - (67) 

2 I Bt(r2)(T3 At{r2,N2)l2 / 



with 



At{Tk,Nk) - Ao e-r.-t + (i _ g-r.t) (i + 27Vfe) , 
BtiTk) = Boe 



(68) 



Finally, if we assume Fi = r2 = F and Ni — N2 ~ N, then 
the covariance matrix (167 > becomes formally identical to (127 1 
and the corresponding Wigner function reads 



Wt(a,/3) 
with 



exp{-2At(|a|' 



')+2BtiaP + a*f3*)} 



At = 



4^Dct[at] 



4^Dct[(Tt 



(69) 



(70) 



If the IPS process is performed on a TWB evolved in a noisy 
environment with both the channels having the same damp- 
ing rate and thermal noise, then the Wigner function of the 
state arriving at the beam splitters is now given by Eq. ( I69l l. 
and the output state is still described by Eq. ( I54t . but with the 
following substitutions 



Aq ^ At , Bq Bt 



crt ■ 



(71) 



V. CONTINUOUS VARIABLE TELEPORTATION 

The scheme of continuous variable (CV) teleportation is 
sketched in Fig.Q A bipartite state Ws is shared between two 
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FIG. 7: Scheme of the CV teleportation. One of the two modes 
of a shared bipartite state Qs is mixed with the input state Qin at a 
balanced BS, and then a double homodyne detection is performed 
on the two output modes measuring the complex outcome 5. The 
teleported state Qout is obtained displacing by the same amount ^ 
the remaining mode of the shared state and averaging over all the 
possible outcomes. 



parties: one mode of the state is mixed at a balanced beam 
splitter (BS) with the state to be teleported. Win, then double- 
homodyne measurement is performed on the two emerging 
modes. The complex outcome ^ of the measurement is used 
in order to displace the remaining mode of Wg and the tele- 
ported state Wout is obtained averaging over all the possible 
outcomes. Here we address the teleportation of the coherent 
state I a), whose Wigner function reads 



Win{z)^-exp{-2\z-a\^}. 

TT 



(72) 



If we consider the following generic shared state: 

cxp{-F|w|2 - G\w\^ + H{vw + v*w*)} 



Wsiv,w) = 



Tr^{FG~-H 



2^-1 



(73) 

and since the POVM describing the double homodyne detec- 
tion is 



1 



(74) 



(5^^-' {() being the complex Dirac's delta function, the output 
state Wout is given by |2] 



Wout(u)) =7r" / / d'zd'vWiniz) 

X Ws{v,w-OW^{z,v) (75) 
■expl-^^i^l (76) 



TTCTout 



fout 



where 



_ 1 F + G + 2H 



(77) 



in turn, the average fidelity of teleportation of coherent states 
reads as follows: 



(78) 
(79) 



F = 7r / d'wWMWontiw) 
FG-H^ 



FG- + F + G-2H 1 + 2 dout 
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FIG. 8: Plots of the teleportation fidelity "Fjwb assisted by TWB in 
the ideal case (Tt = N — 0, dot-dashed line) as a function of A = 
tanhr. The solid lines are Ftv/b with Tt — 0.2 and, from top to 
bottom, N = 0, 0.1, 0.2, and 0.5. 
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FIG. 10: Plots of the teleportation fidelity Fips assisted by IPS as 
a function of A = tanhr with T = 0.9, e = 1, Ft = 0.1 and 
different values of A'^ (solid lines): from top to bottom N — 0, 0.1, 
0.2, and 0.3. The dot-dashed line is F^^s with T = 0.9, e = 1, and 

rt^N = 0. 




FIG. 9: Plots of the teleportation fidelity Fips assisted by IPS in the 
ideal case (Tt = = 0) as a function of A = tanhr. The dashed 
line is Ftwb, whereas the solid lines are Fips with e = 1 and, from 
top to bottom, T = 1.0, 0.9, 0.8, and 0.5. 



When the shared state is the TWB of Eq. i25i . the average 
fidelity is obtained from Eq. ( I78> with F = G = 2Ao and 
H = 2Bo, i.e., 



i^™B(A) = i(l + A) 



(80) 



whereas in the presence of noise one should use the substitu- 
tions (I71> . F-pfjB is plotted in Fig.|8] When the teleportation is 
assisted by IPS, then the fidelity reads as follows: 



^^IPS — 



1 



E 



Pii{r, T, e) FkGk - Hi + Fk + Gk- 2Hk 



(81) 

withFfe = b-fh, Gk = b-gh, and Hk = 2BoT + hk and all 
the involved quantities are the same as in Eq. i54\ . The results 
are presented in Fig.|51for e = 1 andTt = N = 0. The IPS 
state improves the average fidelity of quantum teleportation 
when A is below a certain threshold, which depends on T (and 
e). Notice that, for r < 0.5, i^[ps(A) is always below i^TOB(A), 
at least for e — 1. The effect of dissipation and thermal noise 
is shown in Fig. 1101 

In order to quantify the improvement and to study its de- 
pendence on T and e, we define the following "relative im- 




FIG. 1 1 : Plots of 7?.F as a function of T and e with r = 0.3; we set 
Tt = N — (upper surface) and Ft = N = 0.1 (lower surface). 
TZf, and, in turn, i^ips, is mainly affected by the IPS transmissivity 
T. 



provement : 



7^J.(r,T,e,^,7V) 



F,ps(r, r,£,r,iv)-F™B(r,r,7V) 



F^,{r,r,N) 



(82) 

which is plotted in Fig.[n] we can see that TZp and, in turn, 
i^ips are mainly affected by T when Ft and N are fixed. In 
Fig- El we plot TZp as a function A = tanhr and the quantity 
7?.lj."^^ defined as follows: 



nP{r,T,e,r,N) 



fiPs(r,r,£j-,jV)-F„e(r,0,0) 
F™B(r,0,Q) 



(83) 

i.e., the relative improvement of the fidelity using IPS in the 
presence of losses and thermal noise with respect to the fi- 
delity using the TWB in ideal conditions (Tt ~ N ~ 0): we 
can see that, for the particular choice of the parameters, not 
only the fidelity is improved with respect the TWB-based tele- 
portation in the presence of the same dissipation and thermal 
noise (solid line in Fig. I12t . but the results can be also bet- 
ter than the ideal case (dot-dashed line). We can conclude that 
IPS onto TWB degraded by dissipation and noisy environment 
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ator on two modes is defined as 13 

n(a,/3) = Daia)i-lf^Di{a) ® D,m-lf'Dl{p) , 

(84) 

where a, (3 G C, a and b are mode operators and Da{ct) 



exp{c 



i} and Di,{P) are single-mode displacement 



operators. Since the two-mode Wigner function W{a, (3) can 
be expressed as L2J 



W{a,l3) 



■n(a,/3) 



(85) 



FIG. 12: Plot of the relative enhancement TZp as, & function of 
A = tanhr with T = 0.9, e = 1, and Ft = Af = 0.1 (sohd line). 
The dot-dashed line is namely, the relative enhancement of the 

fidelity using the de-Gaussified TWB in noisy environment with re- 
spect to the fidelity using TWB in ideal case (see text for details): for 
a suitable choice of the parameters, the teleportation assisted by IPS 
in the presence of dissipation and thermal noise, can have a fidelity 
larger than the one of TWB-assisted teleportation also when this is 
implemented in ideal conditions (i.e., Tt = N = 0). 
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FIG. 13: Plot of the teleportation fidelity as a function of the average 
number of photons A'^ of the shared state in the case of TWB (dashed 
line) and a photon subtracted TWB (solid line) for T — 0.999, 
e = 1, and in ideal conditions (i.e.. Ft = A'^ = 0. The inset is a 
magnification of the region < A'^ < 2. 



can improve the fidelity of teleportation up to and beyond the 
value achievable using the TWB in ideal conditions. 

Finally, in Fig. [O] we plot the teleportation fidelity as a 
function of the average number of photons N of the shared 
state in the case of TWB and a photon subtracted TWB: we 
can see that for a fixed energy of the shared quantum channel 
the best fidelity is achieved by the TWB state. The same result 
holds in the presence of dissipation and thermal noise. 

In the next Sections we will analyze the nonlocality of the 
IPS state in the presence of noise by means of BeU's inequal- 
ities o. 



VI. NONLOCALITY IN THE PHASE SPACE 

Parity is a dichotomic variable and thus can be used to es- 
tablish Bell-like inequalities 1.29,1 . The displaced parity oper- 



Il{a, 0) being the expectation value of I\-{a, (3), the violation 
of these inequalities is also known as nonlocality in the phase- 
space. The quantity involved in such inequalities can be writ- 
ten as follows 

Sdp = n(ai, /3i) + n(a2, + n(ai, /^s) - n(a2, /^s) , 

(86) 

which, for local theories, satisfies \Bup\ < 2. 

Following Ref. [30], one can choose a particular set of dis- 
placed parity operators, arriving at the following combination 
lEl 

BMJ) - ^{Vj, -Vj) + n(-3x/j, -Vj) 

+ sVj) - u{-3Vj, sVj) , (87) 

which, for the TWB, gives a maximum B^p = 2.32 (for J' = 
1.6 X lO^'^) greater than the value 2.19 obtained in Ref. |30]. 
Notice that, even in the infinite squeezing limit, the violation 
is never maximal, i.e., |Sdp| < 2-^2 |32]. 

In Ref. |31] we studied Eq. (jsTJ for both the TWB and the 
IPS state in an ideal scenario, namely in the absence of dis- 
sipation and noise; we showed that, using IPS, the maximum 
violation is achieved for T, e — > 1 and for values of r smaller 
than for the TWB. 

Now, by means of the Eq. (jS^ and the substitutions ( 17 11 1, 
we can study how noise affects Sdp- The results are showed 
in Fig. ^5 for e — 1: as one may expect, the overall effect of 
noise is to reduce the violation of the Bell's inequality. When 
dissipation alone is present (N — 0), the maximum of vio- 
lation is achieved using the IPS for values of r smaller than 
for the TWB, as in the ideal case. On the other hand, one 
can see that the presence of thermal noise mainly affects the 
IPS results. In fact, for Tt = 0.01 and N = 0.2, one has 
I'^d'p^^^I f*^'" ^ range of r values, whereas jSpp^'*! falls 
below the threshold for violation. Note that the maximum of 
violation, both for the TWB and the IPS state, depends on the 
squeezing parameter r. 

In Fig.llSlwe plot 6[)p^' as a function of T and e. We can 
see that the main contribution to the Bell parameter is due to 
the transmissivity T. Moreover, as T ^ 1, the Bell parameter 
is actually independent on e. Note that the values of J' and 
r, which maximize the violation, depend on Ft and N, as one 
can see from Fig.^^ in Fig. we have chosen to fix the 
environmental parameters in order to compare the two plots. 
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even if best results can be obtained maximizing Bj^p 
respect J' and T. 



with 
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FIG. 14: Plots of the Bell parameters Sdp for the TWB (top) and 
IPS (bottom); we set S — 1.6 x 10^^, which maximizes B'^^^\ 
and put T = 0.9999 and e = 1 for the IPS. The dashed lines refer to 
the absence of noise (Tt = N = 0), whereas, for both the plot, the 
solid lines are Sdp with Ft = 0.01 and, from top to bottom. A*' = 
0,0.05, 0.1, and 0.2. In the ideal case the maxima are B- 




2.32 and B 



(IPS) 



(TWB) 
DP 



2.43, respectively. 



VII. NONLOCALITY AND HOMODYNE DETECTION 

In principle there are two approaches how to test the Bell's 
inequalities for bipartite state: either one can employ some 
test for continuous variable systems, such as that described in 
Sec. lVll or one can convert the problem to Bell's inequalities 
tests on two qubits by mapping the two modes into two-qubit 
systems. In this and the following Section we will consider 
this latter case. 

The Wigner function W,ps(q;,/?) given in Eq. ( I54t is no 
longer positive-definite and thus it can be used to test the vi- 
olation of Bell's inequalities by means of homodyne detec- 
tion, i.e., measuring the quadratures and x^p of the two IPS 
modes a and b, respectively, as proposed in Refs. 1 16, 17]. In 
this case, one can dichotomize the measured quadratures as- 
suming as outcome +1 when x > 0, and —1 otherwise. The 
nonlocality of W,ps(q;, /?) in ideal conditions has been studied 
in Ref. 1 3 1 ] where we also discussed the effect of the homo- 
dyne detection efficiency jju . 

Let us now we focus our attention on W,ps(q!, (3) when the 
IPS process is applied to the TWB evolved through the noisy 
channel, namely, using the substitutions illi . After the di- 
chotomization of the homodyne outputs, one obtains the fol- 
lowing Bell parameter 

Bud - E{^i,^i) + E{,h,V2) + E[^2,Vi) - £{^2,^2) , 

(88) 

where i?fe and (pk are the phases of the two homodyne mea- 
surements at the modes a and b, respectively, and 



(IPS) 




FIG. 15: The surfaces are plots of the Bell parameters Bdp for the 
IPS state as a function of T and e for different values of Tt and 
iV = 0: (top) Ft = 0; (bottom) Ft = 0.005. We set J = 1.6 x 10"^ 
and r = 1.16 The value of the Bell parameter is mainly affected by 
T. 



dxi)^ dx^^ sign[a;^^ x^J P(xtf^,x^J 



(89) 

P{x^^ , x^^ ) being the joint probability of obtaining the two 
outcomes x^^ and x^p^ WA- As usual, violation of Bell's in- 
equality is achieved when |Shd| > 2. 

In Fig. we plot Shd for §1 = 0, ^2 = 7r/2, ipi — 
— 7r/4 and Lp2 — 7r/4: as for the ideal case 1171 13111 . the Bell's 
inequality is violated for a suitable choice of the squeezing 
parameter r. Obviously, the presence of noise reduces the 
violation, but we can see that the effect of thermal noise is not 
so large as in the case of the displaced parity test addressed in 
Sec. lVIK see Fig. ll4l l. In Fig.^]we plot Shd as a function of 
T and e: as for the displaced parity test (see Fig.llSll. we can 
see that the main contribution to the Bell parameter is due to 
the transmissivity T. 

Notice that the high efficiencies of this kind of detectors 
allow a loophole-free test of hidden variable theories fsS*], 
though the violations obtained are quite small. This is due to 
the intrinsic information loss of the binning process, which is 
used to convert the continuous homodyne data in dichotomic 
results 



We conclude that, considering the displaced parity test in 
the presence of noise, the IPS is quite robust if the thermal 
noise is below a threshold value (depending on the environ- 
mental parameters) and for small values of the TWB parame- 
ter r. 



Vm. NONLOCALITY AND PSEUDOSPIN TEST 

Another way to map a two-mode continuous variable sys- 
tem into a two-qubit system is by means of the pseudospin 
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FIG. 16: Plots of the Bell parameter Bhd for the IPS states for two 
different values of the homodyne detection efficiency: rjn ~ 1 (top), 
and rjH = 0.9 (bottom). We set e = 1 and T = 0.99. The dashed 
lines refer to the absence of noise (Ft = N = 0), whereas, for both 
the plots, the solid lines are Bud with Ft = 0.05 and, from top to 
bottom, N = 0, 0.05, 0.1 and 0.2. 



function 

£;(a,b) = ((a.S) (b-S)), 
where a and b are unit vectors such that 



(90) 



Si-S = cosdaS,+ sin?9„ {e"^-S- + e~'-^-S+) , (91) 
b • S = cosi?h 5"^ + sindb (e"^'S- + e-'^^oS+j , (92) 

with S± — ^{Sx ± iSy). In the following, without loss of 
generality, we set (pk. — 0. Finally, the Bell parameter reads 

Bps = ^(ai,bi) + S(ai, bj) + ^(aa, bi) - ^(a2,b2) , 

(93) 

corresponding to the CHSH Bell's inequality \Bps\ < 2. In 
order to study Eq. ( I93> we should choose a specific represen- 
tation of the pseudospin operators; note that, as pointed out in 
Refs. 1 37, 38], the violation of Bell inequalities for continuous 
variable systems depends, besides on the orientational param- 
eters, on the chosen representation, since different Sk leads 
to different expectation values of Bps- Here we consider the 
pseu dospin operators corresponding to the Wigner functions 

im 

W,ia) = isign[3fie[a]] , W,{a) = -^5^^\a) , (94) 



TT 

Wyia) 



1 

2^ 



5{dic[a])V- 



1 



(95) 



where V denotes the Cauchy's principal value. Thanks to (|94| 
one obtains 



Bud 




FIG. 17: The surfaces are plots of the Bell parameters Shd for the 
IPS state as a function of T and e for A'^ = and different values of 
Ft: (top) Ft = 0, and (bottom) Ft = 0.025. We set r = 0.82 and 
Vn = 0.9 



test: this consists in measuring three single-mode Hermitian 
operator Sk satisfying the Pauli matrix algebra [Sh, Sk] — 
2iehkiSi, si — 1, h,k,l = 1,2,3, and Shk i is the totally 
antisymmetric tensor with £123 = +1 13611 . For the sake 
of clarity, we will refer to 5*1, ^2 and 5*3 as Sx, Sy and Sz, re- 
spectively. In this way one can write the following correlation 



£^TWB (a, b) = cos da COS -dh 
2 sin sin -dh 



for the TWB, and, for the IPS, 



£;iPs(a,b) 



4 

E 

k= 



Ck 



\ Pii{r, T,e) 
2 sin -da sin -db 



arctan [sinh(2r)] , (96) 



cos da COS db 



TrAk 



arctan 



f2BoT + hk 



(97) 



where Ak = [b - fk){b - gk) - {2BoT + hkf, and all the 
other quantities have been defined in Sec.lIIII 

In Fig. [m we plot Bps for the TWB and IPS in the ideal 
case, namely in the absence of dissipation and thermal noise. 
For all the Figures we set -da^ = 0, ~ "^12, and — 
—db2 = 7r/4. As usual the IPS leads to better results for small 



values of r. Whereas B 



(TWB) 
PS 



2V2 



as r 



b: 



(IPS) 
PS 



has a maximum and, then, falls below the threshold 2 as r 
increases. It is interesting to note that there is a region of small 



values of r for which B 



(TWB) 
PS 



< 2 < B 



(IPS) 
PS 



i.e., the IPS 



process can increases the nonlocal properties of a TWB which 
does not violates the Bell's inequality for the pseudospin test, 
in such a way that the resulting state violates it. This fact is 
also present in the case of the displaced parity test described in 
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FIG. 18: Plots of the Bell parameter Bpg in ideal case (Ft = iV = 
0): the dashed line refers to the TWB, whereas the solid lines refer to 
the IPS with e = 1 and, from top to bottom, T = 0.9999, 0.99, 0.9, 
and 0.8. There is a threshold value for r below which IPS gives a 
higher violation than TWB. Note that there is also a region of small 
values of r for which the IPS state violates the Bell's inequality while 
the TWB does not. The dash dotted line is the maximal violation 
value 2^2. 
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FIG. 19: Plots of the Bell parameter Bps for Ft = 0.01: the dashed 
line refers to the TWB, whereas the solid lines refer to the IPS with 
£ = 1 and, from top to bottom, T = 0.9999, 0.99, 0.9, and 0.8. The 
same comments as in Fig. llSl still hold. 



Sec. I VII but using the pseudospin test the effect is enhanced. 
Notice that the maximum violations for the IPS occur for a 
range of values r experimentally achievable. 

InFig.[l9l we consider the presence of the dissipation alone 
and vary T. We can see that IPS is effective also when the ef- 
fective transmissivity T is not very high. We take into account 
the effect of dissipation and thermal noise in Figs.|20| andl21l 
we can conclude that IPS is quite robust with respect to this 
sources of noise and, moreover, one can think of employing 
IPS as a useful resource in order to reduce the effect of noise. 
In Fig.l22lwe plot B^^l^^ as a function of T and e: the main 
effect on the Bell parameter is due to the transmissivity T, as 
in the precious cases. 



IX. NONLOCALITY AND ON/OFF PHOTODETECTION 

The nonlocality test we are going to analyze is schemat- 
ically depicted in Fig. |23] two modes of the de-Gaussified 
TWB radiation field, a and 6, described by the density ma- 



FIG. 20: Plots of the Bell parameter Sps for different values of Vt 
and in the absence of thermal noise (TV = 0): the dashed lines refer 
to the TWB, whereas the solid ones refer to the IPS with T = 0.9999 
and e = 1; for both the TWB and IPS we set, from top to bottom, 
Vt = 0, 0.01, 0.05, and 0.1. The dash dotted line is the maximal 
violation value 2\/2. 
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FIG. 2 1 : Plots of the Bell parameter Bps for Yt = 0.01 and different 
values TV = 0: the dashed lines refer to the TWB, whereas the solid 
ones refer to the IPS with T = 0.9999 and e = 1; for both the TWB 
and IPS we set, from top to bottom, TV = 0, 0.01, 0.1, and 0.2. 
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FIG. 22: The surfaces are plots of the Bell parameters Spg for the 
IPS state as a function of T and e for TV = and different values of 
Vt: (top) Vt = 0, and (bottom) Vt = 0.025. We set r = 0.86. 



trix Q, are locally displaced by an amount a and /3 respec- 
tively and, finally, they are revealed by on/off photodetectors, 
i.e., detectors which have no output when no photon is de- 
tected and a fixed output when one or more photons are de- 
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FIG. 23: Scheme of the nonlocality test based on displaced on/off 
photodetection: the two modes a and & of a bipartite state g are lo- 
cally displaced by an amount a and (3 respectively, and then revealed 
through on/off photodetection. The corresponding correlation func- 
tion violates Bell's inequalities for dichotomic measurements for a 
suitable choice of the parameters a and /3, depending on the kind 
of state under investigation. The violation holds also for non-unit 
quantum efficiency and non-zero dark counts. 



tected. The action of an on/off detector is described by the fol- 
lowing two-value positive operator-valued measure (POVM) 

{no,^,i3,ni,^,B} 11 



n 



0,ri,D 



k=0 ^ 



D 



\k){k\, (98a) 



(98b) 



7] being the quantum efficiency and D the mean number of 
dark counts, i.e., of clicks with vacuum input. In writing 
Eq. ( I98I I we have considered a thermal background as the ori- 
gin of dark counts. An analogous expression may be written 
for a Poissonian background |19]. For small values of the 
mean number D of dark counts (as it generally happens at 
optical frequencies) the two kinds of background are indistin- 
guishable. 

Overall, taking into account the displacement, the measure- 
ment on both modes a and b is described by the POVM (we 
are assuming the same quantum efficiency and dark counts for 
both the photodetectors) 



n|rH«,/3) = n|:''^)(a)®ni'''^)(/3) , 



(99) 



where h,k ^ 0,1, and U^^'"\z) = D{z)Uh,,^^D {z), 
D{z) — exp |za^ — z*a} being the displacement operator 
and z e C a complex parameter. 

In order to analyze the nonlocality of the state g, we intro- 
duce the following correlation function: 



(100) 



h.k=a 



1 + AIr,,Dia, 13) - 2 [g^,Di(^)+yr,,Dm 



where 



(101) 
(102) 
(103) 



and where (A) = TT[gA] denotes ensemble average on both 
the modes. The so-called Bell parameter is defined by consid- 
ering four different values of the complex displacement pa- 
rameters as follows: 

+ E,,,D{a,P')-E^,D{a',P') (104) 

-Ir,,Di(^'^P')~5n,Dia)-yvMP)} ■ (105) 

Any local theory implies that satisfies the CHSH ver- 

sion of the Bell inequality, i.e., \Brj,D\ < 2ya.a'.d.d' L29l. 
while quantum mechanical description of the same kind of ex- 
periments does not impose this bound. 

Notice that using Eqs. (|98} and ( I101M1031 . we obtain 
the following scaling properties for the functions X^.D(a, (3), 
gn,D{a) and3^,,^£)(/3) 



1 



■Qi-i/(i+D){a) 



(106) 
(107) 
(108) 



where = X,,,o, Gri = Gri,Q, and = J^ij.o- Therefore, it 
will be enough to study the Bell parameter for D — 0, namely 
Bri — B,ffi, and then we can use Eqs. ( ll06t -( fTQ8t to take into 
account the effects of non negligible dark counts. From now 
on we will assume D = Q and suppress the explicit depen- 
dence on D. Notice that using expression ( I105> for the Bell 
parameter the CHSH inequality < 2 can be rewritten 

as 

-1 < X„^z5(a,/3) +X,,z?(a',^)+X„,D(a,/9') 

-X„^z,(a',/3')-e^.D(a)-:V^,D(/3) <0, (109) 

which represents the CH version of the Bell inequality for our 
system |39]. 

In order to simplify the calculations, throughout this Sec- 
tion we will use the Wigner formalism. The Wigner functions 
associated with the elements of the POVM J98> for D = are 
given by I19lll 

P^[no.^](0) = ^ exp{-A„|z|2} , (110) 
M/[ni.^](z) = Wmz) - W[Uo,r,]{z) , (111) 

with = 277/(2 - 77), and T4^[I](z) = tt^^ Then, noticing 
that for any operator O one has 

W[D{a)OD''{a)]{z) = W[0]{z - a) , (112) 

it follows that W[D{a) Eo,^ D^a)]{z) is given by 

W[Dia) Ho,^ D\a)]iz) = W[Uo,^]{z - a) , (113) 
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and therefore 

VF[n^^o'°)(a,/3)](z, w) = W^[no.„](z - a) W[no.^]{w - /?) 

(114) 

(115) 

(116) 



W[Uo,^{a) (E) I] (z, w) = W[Uo,rj]{z ~ a) tt'^ 



Finally, thanks to the trace rule expressed in the phase space 
of two modes, i.e., 

Tr[Oi02]=7r2 / d^z d'^wW[0i]iz,w)W[02]{z,w) , 

(117) 

one can evaluate the functions /9), ^,,(0; ), and 3^,,(/3), 

and in turn the Bell parameter in Eq. (I105> . as a sum of 
Gaussian integrals in the complex plane. 

Let us now consider the TWB i25\ . Since the Wigner func- 
tions of the TWB and of the POVM (|99j are Gaussian, it is 
quite simple to evaluate Tjj{a,f3), Gri{a), and 3^,,(/3) of the 
correlation function (I100> and, then, S^; we have 



1]'^ y/Det[cro] 
exp{ -F^dap 



e,(«)=3^,(a) 



X exp ■ 



(2^Det[cro; 



(a/3 - 



A, 



(118) 



-Bl) 
2A„ 



2{Al 



AoAn 



with 



Hr, 



(2Ao + A„)A^^(r) 



2BoM„ir) 



Mr,ir) 



A2 



A{Al-Bl)+4AoA^ 



A2 



(119) 



(120) 
(121) 

(122) 



In order to study Eq. ( I105> . we consider the parametrization 
a — —[3 = J and a' — —(3' = —^/TlJ (more details are 
given in |3). The parametrization was chosen after a semi- 
analytical analysis and maximizes the violation of the Bell's 
inequality (for 77 — 1). In Fig.|25we plot S,, for 77 = 1: as one 
can see the inequaUty \B^\ < 2 is violated for a wide range 
of parameters, and the maximum violation (Btj = 2.45) is 
achieved when — 0.16 and r = 0.74. The effect of non-unit 
efficiency in the detection stage is to reduce the the violation; 
this is shown in Fig.|25j where we plot as a function of J' 
with r ~ 0.74 for different values of the quantum efficiency. 
Note that though the violation in the ideal case, i.e., 77 = 1, 
is smaller than for the Bell states, the TWBs are more robust 
when one takes into account non-unit quantum efficiency. 

In the case of the state ( I54K the correlation function ( I100> 
reads (for the sake of simplicity we do not write explicitly the 




FIG. 24: Plot of Br, for a TWB as a function of J and the TWB 
squeezing parameter r in the case of ideal (i.e., rj — 1) on/off pho- 



todetection. The maximum violation is Br, 
tained when J = 0.16 and r = 0.74. 



2.45, which is ob- 




0.05 0.1 0.15 0.2 0.25 0.3 



J 



FIG. 25: Plot of Br, for a TWB as a function of J with r = 0.74 
for different values of rj: from top to bottom rj — 1.0, 0.9, 0.85, and 
0.80. 



dependence on r, T and e) 
1 



E,,ia,f3) = 1 + 



pii{r, T,e) 



^Cfe{4XW(a,/3) 



k=l 



2[e(^Ha)+3^W(/3)]}, (123) 



where 



X exp { - G(^-) lap - |/3p + iJ^'^) (a/3 + a*/3*)} , 

(124) 



[Gk [Fu + A 
X exp 



(FkGk - Hi) An I ,2 



Gk {Fk + A^) - Hi 
A„ 



(125) 



[Fk {Gk + A„) - Hi] 7? 



X exp 



{FkGk ~Hl)An I |2 
Fk{Gk + A^)~Hl 



(126) 



with F,)"' = F^,"' (r, T, e), G^''^ = G'^' (r, T, e), and iJ^' 



r(fe) 
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FIG. 26: Plot of S,, for the IPS state with T = 0.9999 and e = 1 as 
a function of J and the TWB squeezing parameter r in the case of 
ideal (i.e., r\ = \) on/off photodetection. The maximum violation is 
Br, = 2.53, which is obtained when ^ = 0.16 and r = 0.39. 



H);\r,T,e) given by 

) = A„ - (Ffe + A,) (r, T, e) , (127) 

= A, - (Gfe + A,) AlW(r,r,e) , (128) 

Uf^ =iikMf{T,T,e), (129) 

A^ 

7\/l(''^frTe) = ^ (130) 

where Fk^b- fh, Gk = b- gn, and i/^ = 2BqT - hk and 
all the involved quantities are the same as in Eq. 

In order to study Eq. ( I105> . we consider the parametrization 
a ~ —(3 = J and a' = —[3' = —^/llj. This parametriza- 
tion was chosen after a semi-analytical analysis and maxi- 
mizes the violation of the Bell's inequality (for 77 = 1) L19i] . 
The results are showed in Figs.|26land|^for rj — 1 and e = 1: 
we can see that the IPS enhances the violation of the inequal- 
ity < 2 for small values of r (see also Refs. 1 10, 15, 31|). 
Moreover, as one may expect, the maximum of violation is 
achieved as T ^ 1, whereas decreasing the effective trans- 
mission of the IPS process, one has that the inequality be- 
comes satisfied for all the values of r, as we can see in Fig. 1271 
for T = 0.6. 

In Fig.|28]we plot S,, for the IPS with T ^ 0.9999, e = 1 
and different rj. As for the TWB, we can have violation of the 
Bell's inequality also for detection efficiencies near to 80%. 
As for the Bell states and the TWB, a 77- and r-dependent 
choice of the parameters in Eq. ( 1 1051 1 can improve this result. 

The effect on a non-unit e is studied in Fig. |29] where we 
plot S,, as a function of T and e and fixed values of the other 
involved parameters. We can see that the main effect on the 
Bell parameter is due to the transmissivity T. 

Finally, the effect of dissipation and thermal noise affecting 
the propagation of the TWB before the IPS process is shown 
inFig.EHl 



V 2.2 
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J 

FIG. 27: Plot of Br, for the IPS state as a function of J with r = 0.39 
for different values of T and e = 1 in the ideal case (i.e., 77 = 1): 
from top to bottom T = 0.9999, 0.90, 0.80, 0.70, and 0.60. 
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FIG. 28: Plot of B,, for the IPS state as a function of J with r = 
0.39, T = 0.9999, e = 1, and for different values of 77: from top to 
bottom r? = 1.0, 0.9, 0.85, and 0.8. 



X. CONCLUSIONS 

We have analyzed in details a photon subtraction scheme to 
de-Gaussify states of radiation and, in particular, to enhance 
nonlocal properties of twin-beams. The scheme is based on 
conditional inconclusive subtraction of photons (IPS), which 
may be achieved by means of linear optical components and 
avalanche on/off photodetectors. The IPS process can be im- 




FIG. 29: Plot of S,, for the IPS state as a function T and e with 
J" = 0.16, r = 0.39, and, from top to bottom, -q = 0.99, and 0.90. 
The main effect on Br, is due to the transmissivity T. 
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FIG. 30: Plot of for the IPS state as a function of J with r = 
0.39, T = 0.99, e = 1, r; = 0.9, and n = 0.1 for different values 
of A^: from top to bottom (solid lines) iV = 0, 0.01, and 0.02. The 
dashed line is B„ with = iV = 0. 



plemented with current technology and, indeed, application to 
single-mode state has been recently realized with high condi- 



[1] S. Braunstein, and P. van Loock, Rev. Mod. Phys. 77, 513 
(2005). 

[2] A. Ferraro, S. Olivares, and M. G. A. Paris, "Gaussian States 
in Quantum Information", Napoli Series on Physics and Astro- 
physics (Bibliopolis, Napoli, 2005); e-print, quant-ph/0503237 

[3] Quantum Interferometry III, Special issue of Fort. Phys. 48 
(2000), F. De Martini, G. Denardo, and L. Hardy Eds.; Quan- 
tum Interferometry IV, Special issue of Fort. Phys. 51 (2003), 
F. De Martini Ed. 

[4] Quantum Communication, Computing and Measurement II, 
P. Kumar, G. M. D'Ariano, and O. Hirota Eds., (Kluver Aca- 
demic, Dordrecht, 2000); Quantum Communication, Comput- 
ing, and Measurements III, P. Tombesi, and O. Hirota Eds., 
(Kluwer/Plenum, Dordrecht, 2001). 

[5] A. Furusawa, et al.. Science 282, 706 (1998); S. L. Braunstein, 
and H. J. Kimble, Phys. Rev. Lett. 80, 869 (1998). 

[6] S. L. Braunstein, and H. J. Kimble, Phys. Rev. A 61, 042302 
(2000); J. Jing et al, Phys. Rev. Lett. 90, 167903 (2003). 

[7] H. Yonezawa, T. Aoki and A. Furusawa, Nature 431, 430 
(2004). 

[8] T. Opatrny, G. Kurizki, and D.-G. Welsch, Phys. Rev. A 61, 
032302 (2000). 

[9] P T. Cochrane, T. C. Ralph, and G. J. Milbum, Phys. Rev. A 
65, 062306 (2002). 
[10] S. Olivares, M. G. A. Paris, and R. Bonifacio, Phys. Rev. A 67, 
032314 (2003). 

[11] 1. Wenger, R. Tualle-Brouri, and P. Grangier, Phys. Rev. Lett. 
92, 153601 (2004). 

[12] M. S. Kim, E. Park, P L. Knight, and H. Jeong, Phys. Rev. A 
71, 013801 (2005). 

[13] S. Olivares, and M. G. A. Paris, 1. Opt. B: Quantum and Semi- 
class. Opt. 7, S616 (2005). 

[14] S. Olivares, and M. G. A. Paris, Phys. Rev. A 70, 032112 
(2004). 

[15] S. Olivares, and M. G. A. Paris, 1. Opt. B: Quantum and Semi- 
class. Opt. 7, S392 (2005) . 

[16] H. Nha, and H. J. Carmichael, Phys. Rev. Lett. 93, 020401 
(2004). 



tional probability fl!']. 

We found that IPS process improves fidelity of coherent 
state teleportation and show, by using several different non- 
locality tests, that it also enhances nonlocal correlations. IPS 
may be profitably used also on nonmaximmally mixed entan- 
gled states, as the ones coming from the evolution of TWB in 
a noisy channel. In addition, the effectiveness of the process 
is not dramatically influenced either by the transmissivity of 
the beam-splitter used to subtract photons, not by the quantum 
efficiency of the detectors used to reveal them. 

We conclude that IPS on TWB is a robust and realistic 
scheme to improve quantum information processing with CV 
radiation states. 
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